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We consider the perturbations of the massive vector field around Schwarzschild, Schwarzschild - de 
Sitter, and Schwarzschild - anti - de Sitter black holes. Equations for a spherically symmetric massive 
vector perturbation can be reduced to a single wave-like equation. We have proved the stability 
against these perturbations and investigated the quasinormal spectrum. The quasinormal behaviour 
for Schwarzschild black hole is quite unexpected: the fundamental mode and higher overtones shows 
totally different dependence on the mass of the field m: as m is increasing, the damping rate of the 
fundamental mode is decreasing, what results in appearing of the infinitely long living modes, while, 
on contrary, damping rate of all higher overtones are increasing, and their real oscillation frequencies 
gradually go to tiny values. Thereby, for all higher overtones, almost non-oscillatory, damping modes 
can exist. In the limit of asymptotically high damping, Reu goes to ln3/(8nM), while imaginary 
part shows equidistant behaviour with spacing Imu) n +i — Imui n = 1/4M. In addition, we have 
found quasinormal spectrum of massive vector field for Schwarzschild-anti-de Sitter black hole. 
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I. INTRODUCTION 

The existence of the charged black hole described by 
the Reissner-Nordstrom solution reflects the fact that a 
black hole can possess a massless monopole vector (elec- 
tromagnetic) field. It gets rid of all higher multipoles 
through radiative processes dominated by quasinormal 
ringing at intermediately late times and by power-low or 
exponential tails at asymptotically late times. 

Since Bekenstein's paper it is well-known that a 
black hole cannot possess even a monopole massive vec- 
tor field. Therefore the black hole has to radiate away 
the massive vector field with some quasinormal frequen- 
cies governing this radiation. Nevertheless, the massive 
vector quasinormal modes of a Schwarzschild black hole 
were not studied so far, and, as we shall show in this 
paper, the problem is qualitivcly different from that for 
a massive scalar field, leading to quite unusual quasinor- 
mal behaviour. First of all, let us briefly review what 
we know about massive scalar and massless vector field 
perturbations. 

The massless vector perturbations of the Schwarzschild 
background was considered for the first time in Q. There 
the Maxwell field perturbations were reduced to a single 
wave-like equation for some gauge invariant function ^ = 
*(r,t), 

Here M is the black hole mass and £ is the multipole 
number. Note, that this equation is valid only for £ > 0, 
while for I — (monopole, or spherically symmetrical 
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perturbations) the Maxwell equations in Schwarzschild 
background do not exhibit dynamical degrees of freedom. 
This signifies about existence of non-radiative electro- 
magnetic monopole hair, i.e. about existence of a black 
hole charge. The quasinormal modes and late time be- 
haviour stipulated by this effective potential were found 
in a lot of papers (see recent papers Pj, H, @, H and 
references therein), and are well-studied. In particu- 
lar, we know that massless vector quasinormal modes 
are qualitatively similar to those of scalar or gravita- 
tional fields 0, except for limit of asymptotically high 
overtones: Reui approaches zero for vector field and is 
ln3 1 8wM for scalar and gravitational fields || . 

On the other hand, the massive term corrects the ef- 
fective potential, and for simplest case of scalar field it 
leads to the wave-equation 

/ 2M\ {£(1 + 1) 2M 2 \ n 

*,r»r-.-*,«- 1 a + ~T + m * = °- 
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Here one can take m — and re-cover the massless case. 
The corresponding quasinormal frequencies were found 
in for massless and in for massive case. Massive 
scalar quasinormal modes proved to show quite peculiar 
properties. Thus when one increases the mass of the 
field m, the damping rates of the QN modes decrease 
strongly, so that existence of infinitely long living modes 
called "quasi-resonances" 01 becomes possible. When 
increasing m, lower overtones, one by one, transform into 
"quasi-resonances" , while all the other higher modes re- 
main "ordinary", i.e. damped (For this to happen 
one needs to deal with relatively large values of m, so 
that in a more realistic picture considering backreaction 
of the scalar particle onto a black hole, existence of such 
quasi-resonances is questionable.). On the other hand, 
in 1992 Coleman, Preskill and Wilczek stated that the 
classical vector monopole field is determined by the mass 
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of the field itself, and not by the mass of the black hole 
[l2|. This stimulated the consideration of the late-time 
behaviour of the monopole massive vector field in the 
Schwarzschild background in |13( where the suggestion 
of [T^ was supported by asymptotic treatment. 

One of the earliest papers dealing with massive vec- 
tor field perturbations was that by Galtsov, Pomerant- 
seva and Chizhov |15| . where it was shown that massive 
particles around black holes have quasi-stationary states 
with hydrogen-like spectrum. That was different from 
behaviour of the Proca field in Coulomb potential, where 
bounded states cannot be formed . In the paper 0] , 
the perturbations equations were deduced for the first 
time, yet, as the system of equations for general value 
of multipole number t cannot be decoupled, the solution 
was obtained in the region far from a black hole 
On the contrary in , the perturbation equations were 
reduced to a single wave-like equation, but only for the 
case of spherically symmetrical perturbations and zero 
cosmological constant. 

We are interested now to know what will happen with 
massive vector perturbations in a black hole background. 
In this case the situation is qualitatively different from 
the known massless vector or massive scalar cases. First, 
the wave equation for monopole massive vector pertur- 
bations cannot be reduced to that one for the massless 
vector field, just because the massless vector field does 
not have radiative monopole. Another distinctive fea- 
ture: the corresponding effective potential is not positive 
definite everywhere outside the black hole, so one must 
check the stability of perturbations. 

The most unexpected feature we have found in the 
present paper is that when increasing the mass of the 
field, the lowest frequency and the higher overtones be- 
haviour is qualitatively different: the fundamental mode 
decreases its damping rate what results in appearance 
of infinitely long living modes, while, on the contrary, 
all higher modes decrease their oscillation frequencies, 
leading to appearance of almost non-oscillatory damping 
modes. 

The paper is organised as follows: in Sec I we deduce 
the wave equation for perturbations of the Proca field in 
the background of Schwarzschild, Schwarzschild-de Sitter 
and Schwarzschild-anti-de Sitter black holes. In Sec. II 
the stability of monopole perturbations is proved. Sec. 
Ill deals with quasinormal spectrum for massive vector 
perturbations of Schwarzschild, and Schwarzschild-anti- 
de Sitter backgrounds, including obtaining of the asymp- 
totically high overtone limit. In the Conclusion we give 
a summary of obtained results. 



II. PERTURBATIONS OF PROCA FIELD IN A 
BLACK HOLE BACKGROUND 

We shall consider here the Schwarzschild black hole so- 
lution with a A - term, i.e. Schwarzschild, Schwarzschild- 
de Sitter and Schwarzschild-anti- de Sitter backgrounds 



in which the massive vector field propagates. The black 
hole metric is given by 

ds 2 = -f{r)dt 2 + fir^dr 2 + r 2 {dd 2 + sin 2 6d<p 2 ), (3) 

where 

„, , / 2M Ar 2 \ 

The vector field is described by a four-potential A^, 
which is supposed to interact with gravitational field min- 
imally, i.e. the field equations are generally-covariant 
analogs of the vector field equations in Minkowskian 
space-time. Therefore, the Proca equations 

Fp - m 2 A» = 0, = - A^ v , (4) 

in curved space-time, read 

^=((4r,„ - Ap^g^^g)^ - m 2 A» = 0. (5) 

From here and on the coordinates t, r, and <j) will be 
designated as 0,1, 2, and 3 respectively. 

With respect to angular coordinates we imply ade- 
quate expansion into spherical harmonics. Then, the field 
perturbations can be described by four scalar functions 
of the radial coordinate and time f em (r,t), h tm (r,t), 
k em (r,t), and a tm (r,t): 

A = .f em (r,t)Y em (6,4), (6) 



A x =h tm {r,t)Y tm {e,cj>), (7) 



A 3 = k £m (r, t)Y tm j(6, 0) - a lm (r, t)sin9Y lmfi . (9) 

Considering eq. (5) with fj, = 0, 1 and substituting 
Eqs. (6)-(9) we arrive at the following equations: 

A(fefr - f m ) - iihf - fr m )r 2 ),rf(r) + m 2 r 2 f m = 0, 

(10) 

A(fc^ m - h tm ) - ({f r m - h^)r 2 ) tfir)- 1 + m 2 r 2 h e,n = 0, 

(11) 

where A = £(£ +1). Here we got rid of the function 
a lm (r, t) , so that the final perturbation dynamic can be 
described by the three independent functions of r and t. 

The other two equations of (5), corresponding to \i = 
2, 3, result in a pair of equations with both even and 
odd spherical harmonics. Let us differentiate (10) with 
respect to r and (11) with respect to t. Then, consider the 
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particular case of spherically symmetrical perturbations. 
Thereby taking I — 0, i.e., implicitly, discarding all terms 
containing derivatives with respect to angular variables, 
and introducing the new function 



B — A rt t — At t r, 

we obtain the following equation: 

B tt (2 2M 4Ar 



(12) 



f(r)B, 



f(r) 



B. 



8M 



2A 

T 



B = 0. 



(13) 



Assuming B ~ e luJt , after introducing of W = Br and 
using of the tortoise coordinate r„: dr* = dr/ f(r), we 
get the radial wave-like equation: 



dri 



w 2 * - V(r)tf = 0, 



with the effective potential 
V(r) = I I 



6M 



(14) 



(15) 



When the A-term vanishes, the wave equation (14, 15) 
reduces to that obtained recently in with the help of 
the Newman - Penrose tetrad formalism. A > 0(A < 0) 
corresponds to asymptotically de-Sitter (anti-de Sitter) 
solutions. 

Yet, for perturbations of general multi-polarity, all four 
equations of (5) can be reduced to the matrix equation 
for three scalar functions x f' Q (r), a — 0, 1, 2, 



dr 2 



(16) 



and the matrix M Q ^(r, uS) cannot be diagonalized by the 
r-independent transformations of the vector ^E' Q ,(r), i.e. 
the set of equations (16) cannot be reduced to the wave- 
like equations (14). 



III. EFFECTIVE POTENTIAL AND STABILITY 

Even though we are limited now by spherically sym- 
metric perturbations, one can hope it is possible to judge 
about stability of the system against massive vector 
field perturbations, because usually, if a system is stable 
against monopole perturbations, it is stable also against 
higher multipole perturbations. The effective potential 
for different values of field mass m is given on Fig. 1., 
2 and 3 for Schwarzschild, Schwarzschild-de Sitter and 
Schwarzschild- anti-de Sitter black holes respectively. To 
prove the stability of perturbations governed by the wave 



equation (14, 15) we need to show, that the correspond- 
ing differential operator 
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d 2 



(17) 



is positive self-adjoint operator in the Hilbert space 
of square integrable functions of r», so that there is 
no normalizable growing solution. This provides that 
all found quasinormal modes are damped. For mass- 
less scalar, vector and gravitational perturbations (as 
well as for a massive scalar perturbations) of a four- 
dimensional Schwarzschild, Schwarzschild-de Sitter, and 
Schwarzschild - anti - de Sitter black holes, the effective 
potential is manifestly positive, and therefore the posi- 
tivity of the self-adjoint operator is evident. As a result, 
the corresponding quasinormal modes for these cases are 
damped. Yet, for massive vector perturbations, as we see 
from Fig. 1., 2 and 3, the effective potential has negative 
values near the event horizon. Nevertheless, the effective 
potential is bounded from below and we can apply here 
the method used in [14|, which consists in extension of 
A to a semi-bounded self-adjoint operator in such a way, 
that the lower bound of the spectrum of the extension 
does not change. For this to perform, let us, following 
|14|. introduce the operator 



d 2 



(18) 



and we know that [l4T |: 



(*,A*) L 2 = -{^*D*) houndary + J dr*(\DV\ 2 + W\y\ 2 ), 

(19) 

where 

W = V+(l- 2 Jf-±pjS'{r)-S\r). (20) 

Thus, we need to find the function S(r) which would 
make the effective potential W positive. After investi- 
gation of the form of the effective potential one can see 
that there is a set of functions S(r) satisfying this re- 
quirement. For instance, the function 



S(r) = i (l - — - ^ 
r V r 3 



creates the following potential: 

4(-3M + r)A - m 2 (6M - 3r - r 3 A) 



W 



3r 



(21) 



(22) 



Using Mathematica, one can show that this potential is 
positive outside the event horizon of a black hole. Thus 
a symmetric operator A is positive definite outside the 
black hole for positive and zero cosmological constant, 
and so is the self-adjoint extension. Yet, for the case 
of asymptotically anti-de Sitter black hole, the range of 
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V(r) 



-0.01 




Figure 1: Effective potential as a function of radial coordinate 
for Schwartzschild black hole for different values of mass of 
the field: m = 0.01 (bottom), m = 0.1, and m = 0.2 (top); 
M = 1. 
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Figure 3: Effective potential as a function of radial coordi- 
nate for Schwartzschild- anti-de Sitter black hole for different 
values of mass of the field: m — 0.01 (bottom), m = 0.1, and 
m = 0.2 (top); M = 1, A = 0.05. 
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Figure 2: Effective potential as a function of radial coordinate 
for Schwartzschild-de Sitter black hole for different values of 
mass of the field: m = 0.01 (bottom), m = 0.1, and m = 0.2 
(top); M = 1, A = 0.05. 



IV. QUASINORMAL MODES 

We shall be restricted here by consideration of quasi- 
normal modes of asymptotically flat and AdS black holes 
as those which are most physically motivated. QNMs of 
asymptotically flat black holes may be observed by future 
generation of gravitational antennas |17|. while asymp- 
totically AdS black holes have direct interpretation in the 
conformal field theory in the regime of strong coupling 

m 

Let us start with asymptotically flat black holes. 
The effective potential (17) approaches constant val- 
ues both at event horizon and spatial infinity in this 
case. Therefore the standard QN boundary conditions 
ij/ ~ e ±ik±r, ^ r * ±oo (k + — u, k- = y/u 2 — to 2 ) are 
reasonable. More accurately, taking into consideration 
the sub-dominant asymptotic term at infinity, the QN 
boundary conditions are 

*(r*) ~ C+e lxr *r (lMm2/x) , {r,r* ->+«>), (23) 



the tortoise coordinate is incomplete. At the same time, 
since the effective potential is divergent at spatial infin- 
ity the Dirichlet boundary conditions \&(r = oo) = is 
physically motivated. Then, the boundary term in (19) 
does not contribute to the spectrum, and we obtain the 
positive self-adjoint extension of A. Thereby, we have 
proved that the Schwarzschild, Schwarzschild-de Sitter 
and Schwarzschild-anti- de Sitter space-times are stable 
against monopole massive vector field perturbations. It 
means that there are no growing quasinormal modes in 
the spectra of these perturbations. In the next section we 
shall compute the quasinormal modes for the asymptoti- 
cally flat and AdS cases, and show, that all found modes 
are damped implying the stability. 



X = V ' LO 2 — TO 2 . 

Note that the sign of x ls to be chosen to remain in the 
same complex plane quadrant as to. 

Following the Leaver method, one can eliminate the 
singular factor from 4', satisfying the in-going wave 
boundary condition at the event horizon and (23) at in- 
finity, and expand the remaining part into the Frobenius 
series that are convergent in the region between the event 
horizon and the infinity (see [Tlj for more details). The 
Frobenius series are: 
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Table I: First ten quasinormal modes for Schwartzschild black 
hole: M = 1, m = 0.01 and m = 0.1 and m = 0.25. 



m = 


0.01 


m = 


0.1 


m = 0.25 


n Re (wo) 


-Im(uo) 


Re (wo) 


-Im(wo) 


Re(wo) -Im(wo) 



0.110523 



0.086079 
0.075725 
0.070401 
0.067068 
0.064737 
0.062991 
0.061619 
3 0.060504 
9 0.059575 



0.104649 
0.348013 
0.601066 
0.853671 
1.105630 
1.357140 
1.608340 
1.859320 
2.110150 
2.360860 



0.121577 
0.082277 
0.074036 
0.069451 
0.066451 
0.064299 
0.062660 
0.061359 
0.060293 
0.059400 



0.079112 
0.344140 
0.599791 
0.853002 
1.105200 
1.356830 
1.608110 
1.859140 
2.110001 
2.360730 



0.222081 
0.062605 
0.065511 
0.064570 
0.063243 
0.062006 
0.060925 
0.059991 
0.059182 
0.058475 



0.012994 
0.325191 
0.592979 
0.849359 
1.102860 
1.355170 
1.606850 
1.858140 
2.109180 
2.360050 




Figure 4: Imaginary part of u) as a function of real part of ui 
for first three overtones for increasing mass: m £ (0.01,0.28) 
for n — (diamond), m G (0.01,0.48) for n — 1 (star), m G 
(0.01, 0.75) for n = 2 (box) 



Substituting (24) into (14) we find a three-term recur- 
rence relation: 

a ai+/3oao = 0, a n a n +i +/3„a n +7„a„_i = 0, n > 0, 

Then, using algebra of continued fractions we can find 
the quasinormal modes as those values of ui for which 



a n -iln 



ttn-27n-l 



Pn-2— Q!ra_37n-2/. 



n+l 



^ + l7^^ + 2 



(26) 



i + 27n + 3/-.. 



As this procedure is described in many papers (see for in- 
stance PJj and references therein), we shall not describe 
it here in details, and, go over directly to the obtained 
results. We shall write u = Reu+ilmu), and the damped 
modes should have Iwjuj < 

First of all, let us look at the table I where first ten 
QNMs are presented for three small (with respect to the 
black hole mass M) values of field mass m. Note, that 



as we consider a vector field minimally interacting with 
gravity, i.e. back reaction of the vector field on the met- 
ric is not considered we cannot consider large values of 
m/M. We see from table I that as the overtone number is 
increasing the difference between QNMs for different val- 
ues of m is decreasing and becomes small even at around 
tenth overtone. Thus, one can conclude (and we check 
this by computing high overtones), that high overtone 
behaviour does not depend on the mass term m coming 
into the effective potential (15), that is in agreement with 
previous study of high overtones for massive scalar field 
in Let us remind, that asymptotic limit of the QN 
spectrum for massive vector field does not reduce to that 
for the massless case, because the effective potential (15), 
does not have physical meaning in the limit m = 0. The 
most unexpected feature of the quasinormal spectrum we 
found (see Fig. 4) is that the fundamental mode shows 
correlation with mass of the field m, totally different from 
all the remaining higher overtones. Thus, as the mass m 
is increasing, the real part of the fundamental mode is 
increasing, while the imaginary part is falling off to tiny 
values, leading thereby to existence of the so-called quasi- 
resonant modes, i.e. of infinitely long living oscillating 
modes ^3 1 E3 • On contrary, the second, third (see Fig. 
4) and higher overtones have their real part decreasing to 
tiny values, and, the imaginary part is growing when the 
mass m is growing. Thus higher overtones can lead to ex- 
istence of almost pure imaginary modes which just damp 
without oscillations. We do not know examples of such a 
different behaviour between the fundamental mode and 
higher overtones, for massless fields of any spin [2(3 or f° r 
massive scalar field 0, [ljj, HD> a t least for asymptoti- 
cally flat or de Sitter black holes. The infinitely long liv- 
ing modes can exist for massive scalar field perturbations 
[Toj , but for all modes 1 n °t only for the fundamental 
one. Note however, that for massless vector perturba- 
tions of asymptotically AdS black holes under Dirich- 
let boundary conditions, the fundamental mode is pure 
imaginary (see for instance [2l| and references therein), 
what represents the hydrodynamic mode in the dual con- 
formal field theory J22|. So, the qualitative difference be- 
tween fundamental and higher overtones is not absolutely 
new phenomena, yet, completely unexpected for asymp- 
totically flat space-times. Note also, that despite the fact 
that we have two tendencies: approaching Reu the con- 
stant value log3/87rM , when n is growing, and at the 
same time approaching zero, when m is growing, there is 
no contradiction: to approach the limit log3/87rM, real 
par t of ui should increase again after some certain n = n c 
|23|. We can observe it on the Table I. (third column), 
where we can observe the "local maximum" of Reuj at 
n = 2. 

The modes in Table I and Fig. 4 were found with the 
help of the above described Frobenius technique. Fol- 
lower overtones, one can use, alternatively, the WKB ap- 
proach suggested in [24j and consequently developed to 
3th [2!| and 6th [2(| WKB orders beyond the eikonal ap- 
proximation. The WKB formula has been used recently 
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Figure 5: Imaginary part of u as a function of real part of ui 
for high overtone behaviour: m = 0.01, M — 1. 



in a lot of papers |27j ] and comparison with accurate nu- 
merical data in many cases [28| shows good accuracy 
of the WKB formula up to the 6th WKB order. Here 
we can compare the results with WKB values, but only 
for the fundamental overtone, because for higher ones: 
n > i = 0, and the WKB method cannot be applied. 

Note also, that an effective potential takes negative 
values near the event horizon, and, the WKB formula 
does not take into consideration "sub" scattering by the 
local minimum of the potential and should not be so 
accurate as in the case of the ordinary positive def- 
inite potential. For example, for m = 0.01 we get 
0.110523 - 0.104649i with the help of the Frobenius 
method, and 0.1195- 0.0871i by WKB formula 0. The 
larger m is, the worse convergence of the WKB method. 
Generally we see that the accuracy of WKB approach is 
not satisfactory here because the WKB formula is actu- 
ally good only for i > n. 

Now let us g o over to asymptotically high overtones. 
It is known |ll|. that the mass term does not change the 
infinitely high overtone asymptotic of the Schwarztschild 
black hole. Thus it is natural to expect that the same 
will take place for a massive vector field. Yet, as there 
are no monopole dynamical degrees of freedom for mass- 
less vector perturbations, we cannot formally take m = 
in the considered effective potential. Therefore, using 
the Nollert's method [29j, we computed numerically high 
overtones for non- vanishing values of m (see Fig. 5). 
From Fig. 5 one can learn that as n is growing, the real 
part approaches ItiZ/8'kM , while the spacing in imagi- 
nary part approaches constant: 



Reuj n 



ln3 
8ttM ' 



Imujr, 



(2ra- 1) 
8M ' 



n -> oo. (27) 



It is different from the asymptotic limit which takes 
place for higher multipole perturbations of massless vec- 
tor field |30|: because for the latter case the real part 
of lu asymptotically approaches zero [3(| . To see better 
that in the obtained plot Reuj approaches ln3/8irM let 
us make fit on values n = 1000, 1500, 2000, ....6000. In 



a similar fashion with Nollert's approach, we see that fit 
in powers of is better then in powers of l/n, and 

gives 



Reu n w 0.04372 



0.04780 0.03353 



oo. (28) 



This is very close to Zn3/(87rM) » 0.04373. When in- 
creasing the number of overtones, the obtained fit is 
closer to Zn3/(87rM). 

Following the arguments of [3(j, it is straightforward 
to reproduce numerically obtained asymptotic (28) in an 
analytical way. For this it is enough to remember that 
the effective potential (15), has the following asymptotic 
behaviour in the origin : 



V(r) 



12M 2 



(29) 



and at the event horizon 

V(r) -> const(r-2M) + 0((r-2M) 2 ), r->2M. (30) 

Therefore the general asymptotic solution near the origin 
is 

^(r*) = ciyW* Ji(u>r») + C2y/oJr\J-i (wr»), r — > 0. 

_ (31) 
Then, repeating all relevant steps of |3(| and taking into 
account that near the event horizon the wave function 
has the following asymptotic 

(r*) ~ e 2M»ln( r .(r)-r»(r=2M)) i f ^ ^ 

and equating the two monodromy (which look similar to 
those in jsJj) one gets (28). 

We see that the high overtone asymptotic (28) is the 
same as for gravitational perturbations. This is easily 
understood, because the effective potential looks like that 
for gravitational perturbations with formally taken I = 1 
plus massive term times f(r). Then, as we have shown 
here for vector and in ^lj for scalar fields, the massive 
term does not contribute in high overtone asymptotic. 

The quasinormal behaviour of SAdS black holes is es- 
sentially dependent on radius of a black hole: one can 
distinguish the three regimes of large (r+ >> R), in- 
termediate (r + ~ R), and small (r + << R) AdS black 
holes. From detailed previous study of massless fields, 
one can learn that QNMs of large AdS black holes are 
proportional to the black hole radius, and therefore to 
the temperature ^^|- QNMs of intermediate AdS black 
holes do not show simple linear dependence on radius 
[l8| . Finally, QNMs of small AdS black holes approach 
normal modes of empty AdS space-time |32j . In the limit 
of asymptotically high damping, QNMs show equidistant 
spectrum with the same spacing between nearby modes 
for different massless fields (scalar, electromagnetic and 
gravitational) plj . 

Using the Horowitz-Hubeny method, we obtain the 
quasinormal frequencies for SAdS black hole numerically. 



7 



Table II: Fundamental quasinormal modes for large (r+ = 100-R) intermediate (r+ = li?), and small (r+ = l/10i?) 
Schwartzschild-anti-de Sitter black hole. 







r+ = WOR 


r+ = 1R 


r+ = 1/10J? 


m 


Re(uj( 


i) -Im(wo) 


Re(a!o) -Im(o>o) 


Re(a;o) -Im(a;o) 



0.01 184.959733 266.38559 2.798314 2.671325 2.6929 0.1010 

0.05 185.109096 266.671681 2.800496 2.674197 2.6949 0.1012 

0.1 185.571219 267.521739 2.807245 2.683084 2.700 0.103 

0.15 186.325972 268.912194 2.818276 2.697634 2.709 0.1035 

0.2 187.352413 270.806932 2.833289 2.717471 2.7247 0.1039 

0.25 188.625045 273.161253 2.851928 2.742088 2.7416 0.1055 



As this method is described in a lot of recent works, 
we shall outline only the key points of it here. The 
Schwarzschild-AdS metric function can be written in the 
form 

/(»■) = 1-7 + ^. (33) 

where R is the anti-de Sitter radius. The corresponding 
effective potential is divergent at infinity and is polyno- 
mial function of r. Therefore, one can expand the wave 
function \I/ near the event horizon in the form: 

oo 

^{x) = J2 a n(x-x+) n , x+ = l/r+. (34) 
Here r+ is the largest of the zeros of the metric function 



f(r). The Dirichlet boundary conditions, we shall use 
here, imply that 

|*(r = oo)|=0. (35) 

Then we need to truncate the sum (34) at some large n = 
N, in order to observe the convergence of the values of 
the root of the equation (35) u) to some true quasinormal 
frequency. 

The fundamental quasinormal frequencies are shown 
in Table II for large, intermediate, and small SAdS black 
holes for different values of m. From Table II one can see 
that both real and imaginary parts of the quasinormal 
frequency are increasing when the mass of the field is 
growing. 



Finally, let us find the normal modes of pure AdS 
space-time for the case of massive vector field. The met- 
ric function f(r) of pure AdS space-time has the form: 



/M-l + jjr 



(36) 



We can put the anti-de Sitter radius to be R = 1 in 
further calculations. The tortoise coordinate is connected 
with the Schwarzschild radial coordinate by the relation: 



r = tanr . 

Then, the effective potential has the form: 

2 to 2 



V 



(37) 



(38) 



Let us introduce a new variable 



• 2 * 

z — sin r . 



(39) 



Then the wave equation can be written in the form: 

2 m' 



4z(l-z)* zz (z)+2(l-2z)* z (z)+(^ 2 - - 
After introducing a new function 

* = 9z a (l-z) fi , 
the wave equation takes the form: 



1 



* = 0. 
(40) 



z(l - z)$, zz {z) + ( - + 2a - (2a + 2/3 + l)z ) $, z (z)+ 
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Table III: 


Higher overtones for larf 


;e (r+ = 100i?) Schwartzschild-anti-de Sitter black hole. 






m 


= 0.1 




m = 0.2 


m = 0.3 
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Re(uj ) 


-Im(wo) 


Re(tJo) 


-Im(o;o) 


Re(a;o) 


-Im(tJo) 


1 


316.780 


492.755 


318.601 


495.990 


321.437 


501.024 


2 


447.102 


717.863 


448.936 


721.089 


451.792 


726.109 


3 


577.201 


942.921 


579.043 


946.141 


581.913 


951.155 


4 


707.218 


1167.952 


709.064 


1171.171 


711.940 


1176.185 


5 


837.194 


1392.973 


839.044 


1396.191 


841.925 


1401.203 


6 


967.150 


1617.987 


969.002 


1621.204 


971.886 


1626.215 


7 


1097.093 


1842.997 


1098.946 


1846.214 


1101.832 


1851.226 


8 


1227.027 


2068.006 


1228.881 


2071.222 


1231.77 


2022.23 


9 


1356.956 


2293.009 


1358.811 


2296.230 


1361.70 


2301.24 



^ ~ + a - 1 - ^ - 1 1. + tl (-V 2 ) -j- — — — (ot -\- pA * = 0. (41) 
2.2- 2(1 Zj 4 J 

i 

We are in position now to choose the values of a and vanish. The general solution is 
(3, so that the terms proportional to 1/z and 1/(1 — z) 



* = C x z^- a {\ z)lFA - « + - ~, - - a + /3 + ^ " 2 «>*)+ 



+C 2 z a (l- z^F^a + p- 



+ (3+^,^ + 2a,z) 



(42) 



where constants C%, C2 may be complex. We require 
regularity of the solution at the origin z = 0, C\ = 0, 
and, the vanishing of the wave function at spatial infinity 
implies 

a Fi(a + /8-| s a + /9 + |,i + 2a,«) = 0. (43) 

Note also, that the choice a = 1, (3 = j(l — 
jVl + 4m 2 ) corresponds to the above boundary condi- 
tions both at infinity and at the event horizon. 

Now, it is not hard to see that oj has the form: 



cu n - ^/\A\/3 \^2n + 3+-(l--VT+ 4m 2 ) J . (44) 

This is different from the AdS normal modes for mass- 
less vector field u> n — -\/A/3(2n + 2 + £) [3fl |. where 



I = 1, 2, i.e. monopole perturbations are not dynami- 
cal. 

Note, that it is expected that similar to scalar field 
behaviour [32l | , the massive vector quasinormal modes of 
SAdS black holes should approach their pure AdS values 
(44) as the mass of the black hole goes to zero. For the 
fundamental mode, we see that according to formula (44), 
for m = 0.1, one has u> = 3.004, and, according to the 
extrapolation of the data in Fig. VI obtained numerically 
with the help of Horowitz - Hubcny method, lo, indeed, 
approaches some constant value close to 3. Unfortunately 
we cannot check the accurate numerical correspondence 
to the formula (44), because the series (34) converges 
very slowly for small black hole radius, and therefore one 
needs enormous computer time to achieve the regime of 
very small black holes. 



Numerical data for high overtones, in the regime of 



large black holes, is shown in Table III. There one can 
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Re (u) + 



2 . 6 
2 . 5 
2 . 4 



Figure 6: Real part of as a function of the black hole radius 
r + for small AdS black holes: rn = 0.1, R = 1. 



see that, indeed, in concordance with analytical formula 
(34), one has 



1.29-2.25i. 



(45) 



At a sufficiently high n, the above formula is valid in- 
dependently of the value of the mass field to. It is also 
valid for any large r + because in the regime of large black 
holes the QNMs are proportional to the black hole ra- 
dius r + for massive fields as well. One could say that a 
quasinormal mode at high damping consist of two con- 
tributions. One is proportional to an overtone number n 
and thereby equals to a spacing between nearby modes; 
it is called "gap" . Another contribution does not depend 
on n in the limit n — > oo, called "offset". Thus, one has 

u„ = [offset] + [gap]n, [gap] = 2\/37ri? 2 e~ ra/3 /9, 

(46) 

where R is the anti-de Sitter radius, i.e. at high over- 
tones, the spectrum is equidistant with spacing which 
does not depend on to , and is the same as for gravita- 
tional or massless vector perturbations. This should be 
true also for intermediate and small AdS balck holes, yet 
to check this numerically one needs considerable com- 
puter time, because of the slow convergence of the series 
for small black holes. 

All numerical computations in this paper were made 
with the help of Mathematica. 



computations of the QNMs spectrum, which is done for 
Schwarzschild and Schwarzschild-AdS black holes. Quite 
unexpected property, we found, is that the behaviour 
of the fundamental mode and all higher overtones (for 
asymptotically flat case) are qualitatively different: when 
increasing the field mass m, the damping rate of the fun- 
damental mode goes gradually to zero, leading to appear- 
ing of infinitely long living mode, while all higher over- 
tones, on contrary, decrease their Reu what results in 
existence of almost pure imaginary modes, i.e. damping 
modes without oscillations. 

Asymptotics of infinitely high overtones for 
Schwarzschild balck holes are the same as for cor- 
responding gravitational (massless) perturbations. In 
particular, for Schwarzschild black hole, real oscillation 
frequency approaches Zn3/87rM, while damping rates 
become equidistant with spacing equal 1/4M. This 
value of high damping asymptotic, which coincides 
with that for massless scalar and gravitational fields for 
Schwarzschild black holes can be easily explained by 
two factors: 1) the mass term does not contribute to the 
limit of infinite damping of the quasinormal spectrum, 
and 2) when formally taking the limit to = in the 
effective potential which governs the evolution of massive 
vector perturbations, one has the potential which looks 
qualitatively like that for gravitational perturbations. 

For asymptotically AdS black holes the quasinormal 
spectrum is equidistant at high overtones with spacing 
which does not depend on the mass of the field. 
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V. CONCLUSION 



Fortunately the monopole perturbations of the Proca 
field in the Schwarzschild-(A)dS black hole background 
can be reduced to the wave-like equation with some ef- 
fective potential. Even though the effective potential is 
not positive definite everywhere outside black hole, we 
have proved that spherically symmetrical perturbations 
of massive vector field is stable, i.e. there are not grow- 
ing modes in this case. This is confirmed by numerical 
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